Effect of zero energy bound states on macroscopic quantum tunneling in high-T c 

superconductor junctions 
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The macroscopic quantum tunneling (MQT) in the current biased high-T c superconductor Joseph- 
son junctions and the effect of the zero energy bound states (ZES) on the MQT are theoretically 
investigated. We obtained the analytical formula of the MQT rate and showed that the presence of 
the ZES at the normal/superconductor interface leads to a strong Ohmic quasiparticle dissipation. 
Therefore, the MQT rate is noticeably inhibited in compared with the c-axis junctions in which the 
ZES are completely absent. 
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Great attention has been attracted to theoretical and 
experimental studies of the effect of the dissipation on 
the macroscopic quantum tunneling (MQT) in supercon- 
ductor Josephson junctions A In current-biased Joseph- 
son junctions, the macroscopic variable is the phase dif- 
ference between two superconductors and measurements 
of the MQT are performed by switching the junction from 
its metastable zero-voltage state to a non-zero voltage 
state (see Fig. 1 (c)). Heretofore, experimental tests of 
the MQT have been focused on s-wave (low-T c ) junc- 
tionsiSi 3 . This fact is due to the naive preconception that 
the existence of the low energy nodal-quasiparticle in the 
d-wave order parameter of a high-T c cuprate supercon- 
ductos 4 ^ may preclude the possibility of observing the 
MQT. 

Recently we have theoretically investigated the effect of 
the nodal-quasiparticle on the MQT in the d-wave c-axis 
junctions^ (e.g., the Bi2212 intrinsic junction^ and the 
cross- whisker junctions^). We have showed that the ef- 
fect of the nodal-quasiparticle gives rise to a super-Ohmic 
dissipation^*^ and the suppression of the MQT due to 
the nodal-quasiparticle is very weak. In fact, recently, 
Inomata et al. have experimentally observed the MQT 
in the Bi2212 intrinsic junctions. 11 They have reported 
that the effect of the nodal-quasiparticle on the MQT 
is negligibly small and the thermal-to-quantum crossover 
temperature is relatively high (~ IK) in compared with 
the case of s-wave junctionsiAii 

In this paper, we will investigate the MQT in the 
d-wave junctions parallel to the afr-plane (see Fig. 1), 
e.g., the YBCO grain boundary junctionsi2*iii and the 
ramp-edge junctions^ In such junctions, the zero energy 
bound states (ZES)i£*i§iiLi£ are formed near the inter- 
face between superconductor and the insulating barrier. 
(Note that in the do/do junction (Fig. 1 (a)) no ZES are 
formed as will be mentioned later.) The ZES are gener- 
ated by the combined effect of multiple Andreev reflec- 



tions and the sign change of the d-wave order parameter 
symmetry, and are bound states for the quasiparticle at 
the Fermi energy. Therefore the ZES may give rise to a 
strong dissipation for the MQT. Note that recently Amin 
and Smirnov have theoretically calculated the decoher- 
ence time of a d-wave qubit and discussed the effect of 
the ZES on the qubit operation^ They, however, phe- 
nomenologically assumed that the system coupled to an 
Ohmic heat bath. Instead, we will calculate the effec- 
tive action starting from microscopic Hamiltonian with- 
out any phenomenological assumptions. Moreover by us- 
ing this effective action we will derive the theoretical for- 
mula of the MQT rate and discuss the influence of the 
ZES on the MQT. 



The grandcanonical Hamiltonian of the d-wave junc- 
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FIG. 1: Schematic drawing of the d-wave Josephson junction, 
(a) do /do junction and (b) do/d^/4 junction, (c) Potential 
U((f>) v.s. the phase difference <f> between two superconduc- 
tors. Uo is the barrier height and ui p is the Josephson plasma 
frequency. 
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tions (Fig. 1 (a) and (b)) is given by 
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where /i is the chemical potential. This Hamiltonian de- 
scribes conduction electrons in a potential W(r) which 
includes a boundary. The second term in H. describes 
the anisotropic attractive interaction of strength g(r). 
This model also includes the insulating tunnel barrier I 
where g(r) = by a suitable choice of W(r). Below, we 
will assume that the tunnel barrier is given by a delta 
function potential, i.e., W(r) — woS(x) and consider the 
high barrier limit zq = mwo/h 2 kp 3> 1 (m is the mass of 
the electron and kp is the Fermi wave number.) which 
corresponds to typical experimental situations. 

The partition function of the system can be written as 
a path integral over the Grasmmann fields^ 
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where (3 = 1/ksT and the Lagrangian is given by 
C[t] = J2a I drtp cr (r,T)d T tpa(r,T) + H[t], In order to 
write the partition function as a functional integral over 
the macroscopic variable (the phase difference (j>), we ap- 
ply the Hubbard-Stratonovich transformation. This in- 
troduces a complex order parameter field A(r,r). Now 
the integrals over the Grassmann fields and |A| can be 
performed by using the Gaussian integral and the saddle 
point approximation, respectively. Moreover, by follow- 
ing the method of Schon and Zaikin, 21 we obtain the par- 
tition function as Z = J T>4>(t) exp (— S e ff[(j)]/h), where 
the effective action S e f / in the high barrier limit is given 
by 
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In this equation, M — C (h/2e) is the mass (C is the 
junction capacitance) and the potential U (0) is described 
by 



d\ <t>Ij{\4>) ~ 4> lext 



(4) 



where Ij is the Josephson current and I ext is the external 
bias current, respectively. The dissipation kernel a(r) is 
related to the quasiparticle current I qp under constant 
bias voltage V by 
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In this equation, -D^(t) is the Matsubara Green's func- 
tion of a free boson 
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where v n — 2m: /hp is the bosonic Matsubara frequency 
(n is an integer). 

Below, we will derive the effective action for the two 
types of the ci-wave junction (do /do and do/d^/^) in or- 
der to investigate the effect of the ZES on the MQT (see 
Fig. 1). In the case of the do /do junction, the node-to- 
node quasiparticle tunneling can contribute to the dissi- 
pative quasiparticle current. However, the ZES are com- 
pletely absent. These behaviors are qualitatively identi- 
cal with that for the c-axis Josephson junctions A On the 
other hand, in the case of the do/d^/^ junction, the ZES 
are formed around the surface of the right superconduc- 
tor dn/4- Therefore the node to the ZES quasiparticle 
tunneling becomes possible. 

Firstly, we will calculate the potential energy U in the 
effective action (3). As mentioned above, U can be de- 
scribed by the Josephson current through the junction. 
In order to obtain the Josephson current we start from 
the Bogoliubov-de Gennes (B-dG) equation^ 
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where tp is the phase of order parameter, u(v) is the 
amplitude of the wave function for the electron (hole)- 
like quasiparticle, h(r) — — h?V 2 /2m — \i + woS(x), and 

A(r - r') = fT 1 J2k A k ex P \- ik ' ( r ~ r ')] is the vol_ 
ume of the superconductor). In the superconductor re- 
gions (do and d n u), the B-dG equation (7) can be trans- 
formed into the eigenequation 
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where, £u = h 2 k 2 /2m + h 2 p 2 /2m, — \i (p = 2nn/D and D 
is the width of the junction). The amplitude of the or- 
der parameter A^ is given by A cos 29 = Ad (9) for d 
and Aosin20 = A<i , 4 (0) for where cos 9 = k/kp. 

The Andreev reflection coefficient for the electron (hole)- 
like quasiparticle r^ e (r e h) is calculated by solving the 
eigenequation (8) together with the appropriate bound- 
ary conditions. By substituting rh e (r e h) into the formula 
of the Josephson current for unconventional superconduc- 



tors (the Tanaka-Kashiwaya formula) 
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we can obtain <f> dependence of the Josephson current. 
Here, A± = A( ±fc;P) , fl± = ^/(huj n ) 2 - |A±| 2 , uj n = 



3 



(2n + l)%/j3H is the fermionic Matsubara frequency. In 
the case of low temperatures -C A ), we get 



I\ sin <fi for 
—I2 sin 20 for 
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where 7 X = 37rA /10ei?Ar, h = w 2 h(3Al/35e 3 N c R% 
(Rn = SttHz 2 1 '2e 2 N c is the normal state resistance of the 
junction and N c is the number of channel at the Fermi 
energy). Note that Eq. (10) is only valid in the case of 
the high barrier limit (zq 1). By using Eq. (4), we 
obtain the analytical expression of the potential U, i.e., 
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As in the case of the 5-wav o 25 i 26 and the c-axis junctions j& 
f7 can be expressed as the tilted washboard potential (see 
Fig. 1(c)). 

Next we will calculate the dissipation kernel a(r) in the 
effective action (3). In the high barrier limit (zo 3> 1), 
the quasiparticle current is given b y 17 i 18 



2g f°° 
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where ijv ~ cos9/zq is the transmission coefficient of the 
barrier I, N^b)(E, 9) is the quasiparticle surface density 
of states (L = d and R = d or d^u), and f(E) is 
the Fermi-Dirac distribution function. The explicit ex- 
pression of the surface density of states is obtained by 
Matsumoto and Shiba. 2 — In the case of do, there are no 
ZES. Therefore the angle 9 dependence of N dg (E,9) is 
the same as the bulk and is given by 



N do (E,9) =Re 
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On the other hand, Nd„ /4i (E, 9) is given by 



Re 



E 2 -A d ^(9) 2 



E 



-7r\A d (9)\S(E)(U) 



The delta function peak at E = corresponds to the 
ZES. Because of the bound state at E = 0, the quasipar- 
ticle current for the do/d^/4 junctions is drastically differ- 
ent from that for the do /do junctions in which no ZES are 
formed. By substituting Eqs. (13) and (14) into Eq. (12), 
we can obtain the analytical expression of the quasipar- 
ticle current in the limit of low bias voltages (eV <C Ao) 
and low temperatures (/? _1 <Ao) as 
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Therefore from Eq. (5), the asymptotic form (r ^> H/Aq) 
of the dissipation kernel for the d-wave junctions is given 
by 



a(r) 
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for do /do 
for do/dir/4 
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where Rq = h/Ae 2 is the resistance quantum. These 
results (Eqs. (15) and (16)) indicate that in the case 
of do /do junctions, the dissipation is the super-Ohmic 
type as in the case of the c-axis junctions^ This can be 
attributed to the effect of the node-to-node quasiparti- 
cle tunneling. Thus the quasiparticle dissipation is very 
weak. On the other hand, in the case of the do/d^/A 
junctions, the node-to-ZES quasiparticle tunneling gives 
the strong Ohmic dissipation which is similar to that in 
normal junctions (a(r) ~ 1/t 2 )i2£*2& Therefore the dissi- 
pation for the do /do junctions is enormously weaker than 
that for the do/d^u junctions. 

Let us move to the calculation of the MQT rate T for 
the d-wave Josephson junctions. At the zero temperature 
T is given by T = lim^^//?) \nZ± Within the WKB 
approximation the partition function Z is evaluated by 
the saddle-point approximation. Then T in the low vis- 
cosity limit is obtained by a perturbative treatment^ 
as r f» Aexp (— Sb/K), where Sb = S e ff[4>B] and <pB 
is the bounce solution. Following the above procedures, 
we finally obtain the central results of this paper, i.e., 
the analytical formulae of the MQT rate for the d-wave 
junctions as2S 
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where c = J™ dy y 4 ln(l + l/y 2 )/ sinh 2 (7ry) w 0.0135, 
C(3) « 1.20 is the Riemann zeta function, x — I e xt/Ii{2), 
and To = 12uj p ^j3Uo/2ithjj p exp (— 36J7o/5fiwp) is the 
MQT rate without the dissipation (Uo is the barrier 
height of the potential U and uj p is the Josephson plasma 
frequency. Note that Uo and u p depend on the type of 
the junction). In Eq. (17) SM is the renormalized mass 
due to the high frequency components (to > uj p ) of the 
quasiparticle dissipation and is given by 
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3 h 2 R Q 
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where K\ is the modified Bessel function. 

In order to compare the influence of the ZES and the 
nodal-quasiparticle on the MQT more clearly, we will 



4 



estimate the MQT rates (17) and (18) numerically. For a 
mesoscopic bicrystal YBCO Josephson junctional (A = 
17.8 meV, C = 20 x 10~ 15 F, R N = 100 QAx = 0.95), 
the MQT rate is estimated as 

_T _ r 83% for d /d 
T ~ 1 25% for d R/4 ' 

As expected, the node-to-ZES quasiparticlc tunneling in 
the do/dir/4 junctions gives strong suppression of the 
MQT rate in compared with the do /do junction cases. 

In conclusion, the MQT in the high-T c superconduc- 
tors has been theoretically investigated. The node-to- 
node quasiparticle tunneling in the d /do junctions gives 
rise to the weak super-Ohmic dissipation as in the case of 
the c-axis junctions^ For the do/d^u junctions, on the 
other hand, we have found that the node-to-ZES quasi- 
particle tunneling leads to the strong Ohmic dissipation. 
We have also analytically obtained the formulae of the 
MQT rate which can be used to analyze experiments. 

In the context of an application to the <i-wave phase 
qubit jlSiSLSSiSSi it is desirable to use the do /do or the c- 



axis Josephson junctions in order to avoid the strong 
Ohmic dissipation. However the effect of the ZES can 
be abated by several mechanisms (e.g., by applying mag- 
netic field or by a disorder in the interface) ^LA£ Therefore 
it is interesting to investigate the MQT of the do/d v /4 
junction in such situations. 

Finally, we would like to comment about a recent ex- 
perimental research. Bauch et al. have succeeded to ob- 
serve the MQT in a YBCO grain boundary bi-epitaxial 
Josephson junction^ In their junction, however, the dis- 
sipation is mainly due to not intrinsic mechanisms (the 
ZES and the nodal-quasiparticles) which discussed in this 
paper but an extrinsic impedance in parallel to the junc- 
tion. Therefore the development of the junction fabrica- 
tion and the measurement techniques will enable us to 
directly compare our theory with experimental results. 
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discussions. This work was partly supported by NEDO 
under the Nanotechnology Program. 
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